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INTERACTION  OF  RELATIVISTIC  ELECTRON  BEAMS 
WITH  FUSION  TARGET  BLOW-OFF  PLASMAS 

The  development  of  high~power  pulse  generators  has  created  interest  in 
the  use  of  relativistic  electron  beams  to  compress  and  heat  small  masses  of 
deuterium  and  tritium  to  fusion  in  a  manner  similar  to  that  proposed  using 
lasers. In  order  to  optimize  thermonuclear  yield,  fusion-pellet  designs 
employ  thin  shells  of  high-atomic-number  material  in  which  the  beam  deposits 
energy. ^5^  Hydrodynamic  modeling  of  pellet  implosion  has  neglected  both  the 
effects  of  the  electromagnetic  field  and  scattering  collisions  in  the  beam- 
heated  plasma  blown  off  of  the  shell. This  work  presents  a  formalism  which 
allows  one  to  determine  the  character  of  beam  deposition  in  the  high-atomic- 
nuiriber  plasma  and  shell  when  these  effects  are  included.  The  Boltzmann 
equation  describing  the  beam  with  a  relativistically-correct  Fokker-Planck 
collision  term  is  solved  with  the  assumption  that  the  elastic-scattering 
collision  time  is  the  shortest  characteristic  time  of  the  system.  (This 
approximation  is  valid  for  cases  of  interest:  the  interaction  of  a  I-3  MeV 
electron  beam  of  about  10  nsec  duration  with  initially-solid  shells  of  heavy 
material.)  One-dimensional  solutions  are  then  calculated  in  two  limiting  cases 
of  interest.  From  these,  the  energy-deposition  profile  and  efficiency  of 
energy  coupling  from  beam  to  plasma  are  determined.  Finally,  the  constraints 
placed  on  beams  for  fusion  in  light  of  the  present  analysis  are  discussed. 

The  equation  describing  the  momentum  distribution  function  of  relativistic 
electrons  interacting  with  a  cold,  high-atomic^number  plasma  may  be  written^ 

where  =  1  +  p^/(mc)^.  The  quantities  u  and  u_  are  energy-dependent 

a  Cl 
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scattering  and  energy-loss  frequencies 
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where 


i 


V/  y  (6) 

2Ug  y  p  Bp  my  o  y 

The  equation  governing  is  obtained  by  setting  secular  terms  in  the  second- 
order  equation  equal  to  zero 


0 


(10) 


Lack  of  particle  conservation  is  due  to  beam  electrons,  slowed  by  dynamic 
friction  to  very  low  energies,  merging  with  the  thermal -electron  backround. 
These  are  not  accounted  for  in  f  since  Eq.  (1)  describes  only  the  distribu¬ 
tion  of  super-thermal  electrons. 

The  rate  at  which  energy  is  transferred  from  the  beam  to  a  unit  volume 
of  plasma  is  obtained  by  taking  the  divergence  of  Eq.  (9).  In  the  steady 
state. 


5- 


(11) 


Q  = 


-V*  q 


eE’§  +  4tt 


00 


Steady-state  solutions  of  /  in  one  dimension  are  now  considered.  The 
plane  x  0  is  chosen  to  divide  a  semi-infinite,  uniform  plasma  occupying  the 
region  x  >  0  from  a  vacuum.  A  monoenergetic,  well-collimated  beam  of  rela¬ 
tivistic  electrons  (particle  flux  propagating  in  the  vacuum  is  normally 
incident  on  the  plasma.  The  plasma  is  assumed  to  be  sufficiently  conductive 
to  exclude  the  vacuum  magnetic  field  associated  with  the  incident  beam  for 
times  of  interest.^  Thus,  a  plasma  return  current  j  =  e$  of  thermal 
electrons  must  flow  towards  the  vacuum-plasma  interface.  The  electric 
field  associated  with  this  current  is 

E  =  Tlj  =  Tle$ 


where  T]  is  the  plasma  resistivity.^  Plasma  electrons  which  reach  the 
interface  are  assumed  to  flow  as  vPXB  surface  currents  along  the  x  =  0 
plane.  This  planar  model  approximates  a  relativistic  electron  beam  inci¬ 
dent  on  a  spherical  shell  of  high-atomic-number  material  for  a  shell 
thickness  small  compared  to  the  radius  of  the  sphere. 

The  beam-electron  distribution  function  in  the  plasma  can  now  be 
written 

^  If 

correct  to  second  order.  The  second-order  term  has  been  simplified  using 
the  secular  equation 


/  =  /qCp.x)  +  |ipA  -  ^ 
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where 


.  _  1  /eE  d _ — \f 

""  2Ug^  p  Bp  mY  BxFo 


and  |Jl  is  the  cosine  of  the  polar  angle  in  a  spherical  momentum-space  geometry 
with  polar  axis  along  x. 

The  beam-electron  flux  is  given  by 


$(x)  =2)7 


When  u  /  and  eEA  are  comparable  in  magnitude,  Eqs.  (12)  -  (l6)  must  be 
E  o 

solved  simultaneously  using  numerical  techniques.  Two  limiting  cases  of 
interest  which  are  amenable  to  analytical  solution  are  now  discussed. 

When  the  plasma  density  is  high  and  the  beam-current  density  is  low, 
dynamic  friction  dominates  over  energy  loss  to  the  electric  field.  Neg¬ 
lecting  electric- field  terms  in  Eq.  (l4)  leads  to  the  solution^*® 


C  X  /  x^ 

Y  ^ 


for  an  incident  monoenergetic  beam  with  y  =  Iti  Eq.  (17) » 


r  ^ 

■/.  *=? 


=  r(Y^)-r(Y) 


r(Y)  = 


SL  =  ,  and  C  is  constant.  The  corresponding  electron  flux  is 
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(20) 


(x)  =  $  exp 


where  =  ^(Yq)  *  r(l),  and  =  $(o).  The  volume  heating  rate  when  E  =  0 
may  be  written 

Qc(x)  =  2? 


f'^O 

I  ®^p(-  2^)  ^ 


The  variation  of  Q_  with  x  is  shown  in  Fig.  1. 

(j 

The  transmission  coefficient  of  beam  current  at  the  interface  can  be 
estimated  by  equating  with  the  positive-going  beam  current  there. 


1  CO  c 


P®/-(P,0)dp  , 

-Q _ +  —  $ 


my 


2  o 


(21) 


(22) 


correct  to  first  order.  With  given  by  Eq.  (17) »  the  transmission 
coefficient  is 


(23) 


For  V  =  3,  T  ranges  from  .75  for  aluminum  to  .to  for  gold. 

^o  C 

When  the  beam  current  is  high  and  the  plasma  density  is  sufficiently 
low,  dynamic  friction  can  be  neglected  in  comparison  to  electric- field 
slowing  down.  In  that  case,  Eq.  (10)  predicts  that  $,  and  therefore  E 
are  constant.  The  solution  of  Eq.  (l4)  to  first  order  with  =  0  for 
an  incident  monoenergetic  beam  is 


/  =  /^  +  ^pA  = 


I  a-p|K(Y*)-K(YQ)- 


6 (y-Y*) 


(24) 
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where  a  and  p  are  constant,  6  =  eE/(mcQ  ), 

o 


and  Y*  (x)=  “  eEx/(mc^).  The  quantity  p  is  determined  by  substitution 

into  Eq.  (l6). 

$  “  ^  m^c^  6  p 

o  3 

The  heating  rate  in  the  absence  of  dynamic  friction  is  given  by  Eq.  (11). 


Qe  =  eE$^ 


Using  Eq,  (22),  the  transmission  coefficient  is 


'^E  (Y^-i)a  +  2ep/3 

O 


A  reasonable  boundary  condition  for  large  x  is  needed  to  determine  a. 
Since  increases  with  x  according  to  Eq.  (2)  and  the  definition  of 
dynamic  friction  must  dominate  over  the  electric  field  for  sufficiently- 
large  x.  The  values  of  and  x  at  which  dynamic  friction  becomes 
inportant  can  be  estimated  by  equating  the  two  terms  on  the  right  side 
of  Eq.  (14) 


3U  /  =  eEA 

e  o 


®  Y*^^  =*=  Y  =  V  -eEx  /(jnc^) 

c  o  c 


=  3€Y^  |a-p[K(Y^)-K(Y^)]| 
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The  region  x  >  x  represents  a  strong  absorber  of  slowed-down  electrons. 
It  is  then  reasonable  to  set  the  negative -going  current  at  x  =  equal 
to  zero.  To  first  order, 

f°  f  1 

$^(x  )  =  2tt  /  pdp</  -  dp  +  —  $  =0 


or 


a  =  p 


I  *  •:(yP-k(y^)  . 

'  r.  J 


(50) 


Simultaneous  solution  of  Eqs.  (29)  and  (50)  yield  and  a.  The  resulting 

variation  of  transmission  with  electric  field  is  shown  in  Fig.  2  for 

V  =5.  The  electric  field  is  related  to  the  penetrating  beam  current  by 
o 


Eq.  (12). 


ef  =  8.4X10"^'^  (Z  +  l)n.9®/^e 

O  ^  ^  L 


(31) 


Here,  0  is  the  plasma  electron  temperature  in  eV. 

The  above  calculations  show  that  although  the  electric  field  reduces 

the  beam  transmission  into  the  plasma,  it  does  increase  the  volume  heating 

rate  above  that  due  to  dynamic  friction  alone.  Its  importance  to*  heating 

can  be  estimated  by  comparing  the  heating  rates  and  averaged  over  x. 

i 

For  y  =3,  Q_/Q  ^  ,  so  that  electric-field  heating  dominates  when 

6  ^  The  region  indicated  by  this  inequality  is  shown  in  Fig.  2.  It  is 

seen  that  electric-field  (or  return-current)  heating  is  important  only  when 
S  is  large  enough  to  reflect  all  but  a  small  portion  of  the  incident  beam. 
Thus,  in  order  to  maximize  energy  transfer  from  beam  to  plasma,  the  plasma 
should  be  prepared  in  a  manner  which  keeps  6  small  everywhere.  In  terms 
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of  the  incident  current  J.  =  e#  /T_,  this  condition  takes  the  form 

1  O  £1 

J.  <  10-1^  02) 

11  cm 

For  incident  currents  of  10®  A/cm^  and  1  keV  temperatures,^  Eq.  (52) 
suggests  that  poor  penetration  occurs  into  plasmas  with  ion  densities 
less  than  10^®  cm”®.  If,  during  pellet  irradiation,  the  blow-off  has 
sufficient  time  to  expand  to  a  low  density  region  of  large  extent,  poor 
coupling  of  the  beam  to  the  dense-plasraa  or  solid  portion  of  the  shell 
is  predicted.  However,  if  the  beam  is  not  strongly  scattered  in  the 
low-density  blow-off,  high  transmission  to,  and  collisional  heating  in 
the  dense  plasma  region  can  occur.  For  blow-off  plasmas  obeying  fluid- 
equation  similarity  solutions,^  the  condition  that  the  beam  not  be 
strongly  scattered  in  the  region  n^^  is 

J.T  <  1O^Y^0/Z^  C/cm^  (55) 

1  o 

where  t  is  the  beam  duration  in  sec.  Relativistic  electron  beams  for 
fusion  should,  and  as  proposed^  do,  satisfy  this  requirement. 

It  should  be  mentioned  that  the  assumption  of  poor  magnetic  field 
penetration  leading  to  the  creation  of  large  electric  fields  in  the 
plasma  is  not  usually  valid  in  present-day,  low- temperature  blow-off 
experiments.  Comparing  the  electromagnetic  skin  depth®  with  the  blow- 
off  thickness^  suggests  that  magnetic  neutralization  occurs  when 

T  >  5X10"®  sec.  (3k) 
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Thus,  temperatures  in  excess  of  about  10  eV  are  required  for  substantial 
return  currents. 

In  summary,  it  has  been  shown  that  the  existence  of  large  return- 
current-generated  electric  fields  in  pellet  blow-off  plasmas  do  not 
increase  beam-energy  deposition  because  they  tend  to  inhibit  beam  penetra¬ 
tion  to  high  plasma  densities. 
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numbers.  Return-current  heating  dominates  to  the  right 
of  the  ^  curve. 


13 


1. 


Director 

Defense  Advanced  Research  Projects  Agency 
Architect  Building 
itoO  Wilson  Boulevard 
Arlington,  VA  22209 

Attn:  Sto.  Ltc.  R,  P,  Sullivan 


2.  Director 

Defense  Nuclear  Agency 
Washington,  D.  C.  20309 

Attn:  DDST,  Mr.  Peter  Haas 
RATN 

APTL,  Technical  Library  (2  copies) 

RAEV 

STVL 

APSI  (Archives) 

3.  Director  of  Defense  Research  and  Engineering 
Washington,  D.C.  20301 

Attn:  DAD  (SK)  Mr.  G.  R.  Barse 

4-.  Commander 

Harry  Diamond  Laboratories 
Washington,  D.C.  20438 

Attn:  AMXDO-RBF,  Mr.  John  Rosaco 
A>IXD0-RBH,  Mr.  S.  Graybill 
AMXDO-RC,  Dr.  Robert  Oswald,  Chief  Lab  300 

5.  Air  Force  Weapons  Laboratory,  AFSC 
Kirkland  AFB,  NM  87II7 

Attn:  DY,  Dr.  Guenther 

EL,  Mr.  John  Darrah 
DYS ,  Dr.  Payton 
SAA 

SUL,  Technical  Library 
ELP,  TRE  Section 

6.  Space  and  Missile  Systems  Organization 
Post  Office  Box  929^0 

Worldway  Postal  Center 
Los  Angeles,  CA  9^009 

Attn:  SKT,  Peter  H.  Stadler 

RSP,  System  Defn.  and  Assessment,  Ltc.  Gilbert 


1 


Sandia  Laboratories 
P.  0,  Box  5800 
Albuquerque,  NM  87 115 

Attn;  Document  Control  for  Dr*  J*  V.  Walker,  5220 
Document  Control  for  52^2,  Dr,  G,  Yonas 
Document  Control  for  Technical  Library 

Aerospace  Corporation 

P.  0,  Box  92957 

Los  Angeles,  CA  9OOO9 

Attn:  J,  Benveniste 

Dr.  Jerry  Comisar 

Austin  Research  Associates,  Inc. 

600  West  28th  Street 
Austin,  Texas  7^705 

Attn;  Dr.  William  E.  Drummond 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  Ohio  45201 

Attn;  P,  Malozzi 

Max^^ell  Laboratories,  Inc. 

9244  Balboa  Avenue 
San  Diego,  CA  92123 

Attn;  Dr.  P.  Korn 

Mission  Research  Corporation 

755  State  Street 

Santa  Barbara,  CA  93101 

Attn;  Dr.  Conrad  L.  Longmire 

Physics  International  Company 
2700  Merced  Street 
San  Leandro,  CA  9^577 

Attn;  Doc.  Control  for  Dr.  Sid  Putnam 
Doc.  Control  for  Mr,  Ian  Smith 

R  &  D  Associates 

P.  0.  Box  3580 

Santa  Monica,  CA  90403 

Attn:  Dr.  Bruce  Hartenbaum 

Science  Applications,  Inc, 

P.  0.  Box  2551 
La  Jolla,  CA  92037 

Attn:  Dr.  J.  Robert  Beyster 


2 


16.  Stanford  Research  Institute 
535  Ravenswood  Avenue 
Menlo  Park,  CA  9^025 

Attn:  Dr.  Robert  A,  Armlstead 

17.  Van  Lint,  Victor  A,  J.  (Consultant) 
7850  Convoy  Court 

San  Diego,  CA  92111 

18.  DDC,  12  cys 

19.  Code  2628,  20  cys 

20.  Code  7700,  25  cys. 

21.  Code  7770,  20  cys 


3 


